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ABSTRACT 

I We investigated the reported distances of Galactic black hole (BH) and neutron star 

low-mass X-ray binaries (LMXBs). Comparing the distances derived for the neutron 
stars Cyg X-2 and XTE J2123-058 using the observed Eddington Hmited photospheric 
radius expansion bursts with the distances derived using the observed radius and 
' effective temperature of the companion star we find that the latter are smaller by 

I approximately a factor of 1.5-2. The latter method is often employed to determine 

. the distance to BH LMXBs. A possible explanation for this discrepancy is that the 

stellar absorption lines in fast rotating companion stars are different from those in 
the slowly rotating template stars as was found before for early type stars. This could 
lead to a systematic mis-classification of the spectral type of the companion star, 
O ' which in turn would yield a systematic error in the distance. Further, we derive a 

j~l ', distance of 4.01^ 2 V404 Cyg, using parameters available in the hterature. 

c/3 ■ The interstellar extinction seems to have been overestimated for XTE J1550-564 and 

possibly for two other BH sources (H 1705-25, and GS 2000+25) as well. As a result 
hereof the distance to XTE J1550-564 may have been underestimated by as much as 
yj ■ a factor three. We find that, using the new distances for XTE J1550-564 and V404 

r> I Cyg, the maximum outburst luminosity for at least 5 but perhaps even 7 of the 15 BH 

^ . soft X-ray transients exceeds the Eddington luminosity for a 10 Mq BH showing that 

these systems would be classified as ultra-luminous X-ray sources had we observed 
them in other Galaxies. This renders support for the idea that many ultra-luminous 
X-ray sources are stellar-mass rather than intermediate-mass BHs. We find that the 
rms value of the distance to the Galactic plane for BHs is consistent with that of 
neutron star LMXBs. This suggests that BHs could also receive a kick-velocity during 
their formation although this has to be investigated in more detail. We find that the 
Galactic neutron star and BH 1 and b distributions are consistent with being the same. 
The neutron star and BH distribution is asymmetric in 1 with an excess of systems 
between -30° < 1 < 0° over systems with 0° < 1 < 30°. 
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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are binary systems in 
which a compact object - either a neutron star or a black 
hole - accretes matter from a low-mass companion star. 
Many of these systems are found to be transient, so called 
soft X-ray transients (SXTs). Several hundre ds have been 
found in our own Galaxy (see lLiu et alJl200lh and many of 
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the X-ray point sources in other Galaxies are likely to be 
LMXBs as well (|^bbian3 [l993) . However, some of these 
X-ray sources in other Galaxies have luminosities in excess 
of the Eddington limits of both ne utron star and t en solar 
mass black hole (BH) LMXBs (cL lFabbianolll989l) . These 
systems are tho ught to contain intermediate mass BHs 
JColbert fc Mushotz kv 1999) or stellar mass BHs where the 
emission is anisotropic ([King et ahlboOlT) or indeed super- 
Eddington jBegelmanll2002ll . 

Determining the distance to LMXBs is important, 
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e.g. to determine the peak or quiescent luminosity and to de- 
termine if the compact object in LMXBs receives an (asym- 
metric) kick velocity at birth or not. Typically, a large dis- 
tance to the Galactic plane implies the occurrence of a kick, 
unless the systems at large z-distances were formed in the 
halo. 

In this Paper we collect distance estimates from the 
literature and discuss possible systematic trends related to 
the used distance estimation methods (Section 2). In Section 
3 we discuss the implications for the Galactic distribution of 
(BH) LMXBs and the Galactic population of ultra-luminous 
X-ray sources (ULXs). 



2 DISTANCES TO LMXBS 
2.1 Black hole distances 

For most BH SXTs it is not feasible to obtain a trigonomet- 
ric parallax measurement. Instead the distance generally is 
determined by comparing the derived absolute V^-band mag- 
nitude with the (dereddened) apparent magnitude, taking 
into account a possible contribution from residual accretion 
(the distance derived using this method is sometimes called 
a photometric parallax). A first guess of the distance can be 
obtained by assuming that the absolute magnitude is that 
of a main-sequence star of the observed spectral type, af- 
ter determining the best-fit spectral type fr om the data, 
e.g. v ia the optimal extraction technique fsee iMarsh et all 
ll994l for a description of this technique). We call this method 
A. However, ideally, the radius, spectral type, and luminos- 
ity class are determined directly from the data. We denote 
this method B. The observations give the orbital period, P, 
the radial velocity of the donor star, K2 and, by comparing 
the donor star spectrum with templates that are Doppler 
broadened, the rotational velocity of the donor (vsin i). The 
inclination can be estimated from modelling the ellipsoidal 
variations in the light curve. From Kepler's third law and 
the assumptions that the donor fills its Roche lobe and is in 
co-rotation with the orbit, all system parameters can be de- 
termined as follows. The mass function and inclination relate 
the masses of the two components, while a combination of 
Kepler's law and the equation for the volume of the Roche- 
lobe gives ^ = QAQ[{l + qfq]^ fe.g. lWade fc HornellQSsI: 
q is defined here as the mass of the secondary divided by the 
mass of the primary), giving an independent estimate of q 
thus providing the system parameters. The radius can be 
estimated from vsini = 2 7rRsim/P. Alternatively, if vsini 
or the inclination are not known, a quite good estimate of 
the radius can be obtained using the fact that the mean den- 
sity of a Roche-lobe filli ng star depends only on the orbital 
period llPaczvhskilll97lH . The radius of the donor star can 
be estimated via R2 = RRochoLobo = 0.234 P^^^ Mj''^ (where 
Porb is in hours, M2 in solar masses, and R2 in units of solar 
radii) . 

In order to estimate the absolute magnitude one either 
uses the surface brightness for the observed spectral type or 
colour , as given e.g. bv .Barnes fc Evans (1976) and Po pper! 
il980l) or uses the determined radius and effective temper- 
ature, together with an appropriate bolometric correction. 
Alternatively, one uses the fact that for given surface bright- 
ness the observed flux scales with the angular diameter (2 



R/d), which together with the radius gives the distance. 
However, there are (small) differences in surface brightness 
for different luminosity classes of the same spectral type. 

In case of GRQ J 1655 -40 and GRS 19 15-H05, 
iHiellming fc RupenI ||1995H and lFender et al.1 lll99gD . respec- 
tively determined limits on the distance from the observed 
proper motion for receding and approaching blobs, assuming 
the jet ejections are intrinsically symmetric and noting that 
the maximum velocity of the ejections is the speed of light. 
We list this method as method C. The distance can also 
be estimated using the interstellar absorption properties of 
the source. There are different ways to do this. For several 
interstellar absorption lines and diffuse interstellar bands 
the observed equ ivalent width correlates with colour excess 
(cf.|Herbi3[l99|). The colour excess can be converted in a 
distance estimate (e.g. using the calibration of lBeals fc Qkd 
Il953h . It is also possible to compare the observed extinction 
to that of (OB) stars in the observed field for which the 
distance is known or can be derived (e.g. Ivan Paradiis et alJ 
[1986). The distance can also be constrained by using high- 
resolution spectroscopic observations of the interstellar ab- 
sorption lines to trace individual gas clouds and their ve- 
locities. Assuming that the apparent velocities of the dif- 
ferent gas clouds are projected Galactic rotation velocities, 
a l ower limit on the distance of the object can be found 
(cf.lH^nesenniiOQi). We refer to distances derived using 
the interstellar absorption properties as method D. Finally, 
we note here that various other methods to estimate the dis- 
tance to BH SXTs have been proposed and used. We do not 
discuss these in det ail but we merely mention some of them: 
iHvnes et al.l (|2002!) used the n ormalisation of a model de- 
scribing the accretion disc fiux, iMaccarond (|200^ used the 
tra nsition between the h igh/soft and low/hard X-ray states, 
and I Jonker et all (l2004ll proposed to use the normalisation 
of the radio - X-ray correlation in the low/hard state as a 
distance indicator. 

In Table Q we show the distance estimates in the lit- 
erature based on method A, B, C, and D for sources for 
which the best estimate of the mass of the compact object 
is above the upper limit o f the neutron star mass of ~ 3 Mq 
iRhoades fc RuffinillTgrl iKaloeera fc Bavml 119961) . Below, 
we briefiy discuss the sources listed in Table^for which dis- 
parate values for the distance exist in the recent literature. 
lA 0620-00: 

[karsh et al. ( 1994^1 report a distance of 485 pc for lA 0620- 
00 whereas most other workers report a distance close to 1.2 
kpc (for references see Table 0. This is due to a numerical 
error of a factor of 2 in the result of Marsh et al. (1994; also 
Marsh private communications). Hence, the distance Marsh 
et al. (1994) find is 912/sini pc. The correct equation for 
the stellar distance is: 

2E ^^ QA6q^/^ jl + qf^'' K2P ^^ Pjvsmi) 
^ e 27r0sini ^ 27r6'sini ^ ' 

Where is the stellar diameter. 
GRQ JQ422+32 : 

IWebb et all J2000ll report a distance of 1.39±0.15 kpc. How- 
ever, besides the fact the best-fit spectral type they derived 
for the companion star was M4- 5, later than that derived 
by e.g. iHarlaftis et"all Jl999l) and lGelino fc HarrisonI (|20o3) 
who found a Ml-2 spectral type, they used method A which 
is known to underestimate the radius of the companion star. 
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Table 1. BH SXT distance estimates. We indicate whether method A, B, C, or D has been used to derive the distance (see text). The 
z-values have been rounded to the nearest 25 pc. 



Name 


1 


b 


Spectral 


^ orb 


Distance Sz Method 


z 


References 








type 


(hours) 


(kpc) 


(pc) 





GRO J0422+32 


165.97 


-11.99 


K9-M2V 


5.09 


2.5-3.0 (B) 


-525 - -625 


[1,2] 


lA 0620-00 


209.96 


-6.54 


K4-K5V 


7.75 


1.2+0.4 (B) 


-125 


[3,4] 


GS 1009-45 


275.88 


9.35 


K6-M0V 


6.84 


5.7+0.7 (B) 


925 


[5,6] 


XTE J1118+480 


157.66 


62.32 


K5-M1V 


4.08 


1.8+0.6 (B) 


1600 


[7,8] 


GS 1124-684 


295.30 


-7.07 


K3-K5V 


10.4 


5.5+1.0 (B) 


-675 


[9,10] 


4U 1543-47 


330.92 


5.43 


A2V 


26.8 


7.5+0.5 (B) 


700 


[11,12] 


XTE J1550-564 


325.88 


-1.83 


G8IV-K4III 


37.0 


5.3+2.3 (S) 


-175 


[13] 


GRO J1655-40 


344.98 


2.46 


F2-F6IV 


62.9 


3.2+0.2 (C) 


125 


[14,15,16,17] 


GX 339-4 


338.94 


-4.33 


?? 


42.1 


>6 (D) 


-450 


[18] 


H 1705-250 


358.59 


9.06 


K0-K5V 


12.5 


8.6+2 (B) 


1350 


[4,19,20] 


SAX J1819.3-2525 


6.77 


-4.79 


B9III 


67.6 


9.6+2.4 (B) 


-800 


[21] 


XTE J1859+226 


54.05 


8.61 


G5-K0V 


9.17 


6.3 + 1.7" 


950 


[22,23] 


GRS 1915+105 


45.37 


-0.22 


K-MIII 


34 (days) 


iiti (c) 


-50 


[24,25,26] 


GS 2000+25 


63.37 


-3.00 


K3-K6V 


8.27 


2.7+0.7 (B) 


-150 


[4,27,28] 


GS 2023+338 


73.12 


-2.09 


G8-K1IV 


6.47 (days) 




-150 


[29,30, this work] 



" A method to estimate the distance which has not been discussed here in detail has been used, see reference 23 for more info. 



R eferences: fll iGelino fc Harrisoi]|2003l [2] iHarlaftis et ai]ll99Sl [3] IShahbaz et al.lll994al [4] iBarret et al.lll996l [5] iBarret et al.ll200(l 
[6] iFilippcnko ct al."l999^. [7] 'McClintock ct al."2001. [81 'Wagner ct al."200l' [91 'Orosz ct al."l996'. [lOl'Esin ct al 199/. [lll'Oros z et alj 
[l99 8. [121 Orosz ct al. 2002a, [13] Orosz ct al. 2002b, [14] Shahbaz ct al. 1999, [15] Orosz & Bailvn 1997, [161 Tingav ct al. 199^ [171 
Tlicllming & Ruocn 1995, [181 Hynes et al . 2004, [191 Romillard ct al. 1996, [201 Harlaftis ct al. 1997, [21] Orosz ct al. 2001, [22] 
iFilippenko fc Chornock.,200L. [23llHvnes et al....2002.. [241 .Mirabel fc Rodrigucz,.1994, [25] .Fender ct al. 1999. [261 Grciner et al..i200U 
[271 IHarlaftis et al.lll996L [281 ICallanan et alj|l996l [291 IWagner et alj|l99a [301 Ishahbaz et al.lll994bl 



In the distance determination lWebb et al] i2000l) use the ab- 
solute maRnitude of a s tar with radius R = 0.24 Rq (from 
iKirkpatrick et aLlllQQ^ via the bolometric magnitude and 
the effective temperature), whereas the radius one would de- 
rive using the equation in the first paragraph of Section [2. II 
would give R = 0.53 R©. 

For that reason, we pre f er the dist ance derived using 
metho d B by Harlaftis ct al. (1999) and lGelino fc Harrison! 
i2003l) . However. iHarlaftis et aL H999I) used the same erro- 
neous e2imtion_j£Marsh_e^ al. (1994). If we correct for this 
error. IHarlaftis et al.' ("1999") found a distance of 5. Oil. 6 kpc. 
Interestingly. iGelino fc Harrison ( 2003i) found a distance of 
2.5±0.3 kpc even t hough they also used method B. This is 
due to the fact that iHarlaftis et al.| l|l999 |) use (V — R) o ~ 1 
for an Ml -2 star, while, both in lEaton fc Poj \l984) and 
ICoxl ('2000) give a value for the {V — R)o colour closer to 1.5 
for a M2V star. Taking a (V - R)o of 1.5 for an M2 as in the 
relation of lEaton fc Pod would yield a distance close 

to 2.5 kpc. Indeed, if we take the radius of the secondary 
and apparent V^-band magnitude corrected for interstellar 
extinction and an ac cretion disc contribution as g iven by 
IHarlaftis et"ai] lll999l) but use the relation given bv lPoppen 
l^8oir "for an M2 star we find a distance of 2.8 kpc. We 
conclude that the distance to GRO J0422+32 is likely to be 
2.5-3.0 kpc. 

GX 339-4: 

In the case of GX 339-4, the distance has been estimated 
to be ~ 4 kpc usi ng the equivalent w idth of the Call-K in- 
terste llar line (e.g. lCowlev et alll987l : see lBuxton fc Vennej 
l2003l for an overview) . Recently, from high resolution spec- 
troscopic observations resolving t he contributions to the in- 
terstellar Na D absorption lines, iHvnes et alj i2004ll found 



that the distance to GX 339-4 is likely to be more than 6 
kpc. They explain that in order for the distance to be ~ 4 
kpc the line-of-sight towards GX 339-4 must be peculiar. 
The limit on the distance of GX 339-4 found by IVIaccaron^ 
(2003) indicates that d>7.6 kpc. In Table [3 we refer to the 
value derived bv lHvnes et al.l ^|2004^ . 

GS 2023+338, a.k.a. V404 Cyg: 

The distance to GS 2023+338 has been reported to 
be close to 3 kpc IShahbaz et alj Il994tfl or 8 kpc 
JWhite fc van Paradii jll996: even a distan ce of 11 kpc has 
been mentioned lCaaE^e^tal ]ll99j.|King (199^ finds that 
if the secondary is a stripped-giant, GS 2023+338 must have 
a distance 3.5-5.1 kpc. In an attempt to reconcile the differ- 
ent distances we recalculate the distance using method B. 
We used the relation between the absolute visual magnitude, 
spectral type and radius of the companion star given by 
IPopperl Jl98Cl) The spectral type of GS 2023+338 is K0±1 
f Casares et al.l 1993). We t ake R2 = 6 ± IRr for t he radius 
of the companion star after IShahbaz et al.l lll994lj) who ob- 
tained this from modelling the ellipsoida l variations in the 
-fsT-band light curve. From the relation of lPoppeJ lll980l) we 
find for the absolute magnitude 1.8 ^ My ^ 3.0. We fur- 
ther take m„ = 18.7 using the observed V^-band magnitude 
of 18.42 and the fact that approximately 25 per cent of the 
light in the F-band was estimated to come from the accre- 
tion disc, althoug h the uncertainty in the accretion disc con- 
tribution is large i Casares et al.lll993l) . Here we took 25±15 
per cent for the accretion disc contribution to the light in the 
F-band; from this we get for the observed 1^-band magni- 
tude 18.5 ^ mv ^ 19.0. Together with an assumed interstel- 
lar absorption of Av ~ 3.3, this yields a distance of 4.01^ 2 
kpc. Note that this calculation does not include an uncer- 
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tainty in the interstellar absorption (see below). Since, as 
mentioned above, the accretion disk contribution is uncer- 
tain we also determined the distance assuming that there 
is no contribution of the accretion disc to the A'-band flux 
(this assumption yields a lower limit to the dist ance). For a 
KOV/III star the {V - K)o colour is 1.96/2.31 l)Co^l20n(ll . 
using the relations between the {V — K)o c olour and the 
surface brightness as determined by iBailevI lll98ll) we get 
a lower limit to the distance of 2.7 kpc for a KOIII and 3 
kpc f or a KOV companion s tar (again we used R2 = 6 ± IRq 
after IShahbaz et al.lir994bl) and we took Ak = 0.4. 

2.2 Neutron star LMXB distances 

In order to compare the BH distances with the neutron star 
LMXB distances, we list in Table |2|the distances to neutron 
star LMXBs. We excluded sources in Globular Clusters since 
we want to compare the neutron star sample with that of 
BHs and (so far) BHs have not been found in Globular Clus- 
ters. 

Some of the type I X-ray bursts, more specifically those 
showing evidence for pho tospheric radius ex pansion can be 
used as a standard candle (Ivan ParadiiJl973) . Using the dis- 
tance of the Galactic Centre and those of Globular Clusters 
to estimate the maximu m peak luminosity f or photospheric 
radius expansion bursts 'van Paradii^ il98ll) found that the 
average peak luminosity of photospheric radi us expansion 
bursts is 3 X 10^* erg s~^. IVerbxmt et al.1 (^3) used a more 
extensive sample of X-ray bursters in Globular Clusters and 
found an aver age peak luminosity of (4.0 ± 0.9) x 10^* erg 
s~^. R.ecentlv. iKuulkers et all (|2003h found a neutron star 
Eddington luminosity of (3.79 ± 0.15) x 10^** erg s"^ for 
the peak luminosity of radius expansion bursts in Globular 
Clusters. This luminosity is consistent with the Eddington 
luminosit y of a 1.4 Mp neutron star accreting helium rich 
material. iKuulkers et al.l lf2003^ also found that a few sys- 
tems have a lower peak luminosity of ~ 2 x 10^* erg s"^ 
which can be interpreted as the Eddington luminosity for 
hydrogen rich accreted material. Therefore, and for reasons 
explained in Section f2. 3. 41 we use both a peak luminosities 
to calculate the distance to the neutron star LMXB. 

For Aql X-1 we determined the burst peak fiux in a 
0.25 s bin from a photospheric radius expansion burst de- 
tected with the 7?Jfri5 satellite in the observation starting on 
MJD 51364.834069 (Terrestrial Time). In this we subtracted 
the average persistent emission 5-105 seconds before the 
burst. The bolometric burst peak flux is 1.1x10"^ erg cm~^ 
s~^. We corrected the Aql X-1 flux for the fact that fluxes 
derived using the RXTE satellite are found to be systemat- 
ically higher by about 20 per cent than the X-ray fluxes in 
the same band found by other satellites (cf. [ Tomsick et all 
Il999bl: iBarret et ai]l2000l: IKuulkers et al.ll2003l) . Finallv. we 
use the distance for Sco X-1 as dete rmined from radio par - 
allax measurements (d=2.8±0.3 kpc. lBradshaw et alll999l) . 

2.3 Systematics and uncertainties in distance 
estimates 

lOrosz fc Kuulker^ (^9^ derived a distance of 7.2±1.1 kpc 
for Cyg X-2, using metho d B (they used the s pectral type of 
the companion derive d bvlCasares et al .119981) . From the ra- 
dius expansion burst JSmal(Jl998ll we flnd a distance of 15.3 



kpc, approximately a factor two larger. The large differences 
in these distance estimates are difficult to explain. Even if 
the radius expansion bursts are hydrogen-rich and hence 
the burst peak luminosity is lower (see Section [2.3.411 . there 
would still be a difference in distance estimates of a factor 
1.5. However, due to the large photon rate deadtime effects 
could lower the apparent burst peak ffux. The bolometric 
burst peak ffux corrected for deadtime effects was 1.5x 10~* 
erg cm~'^ s"'^ (Smale 1998i ; compare this with the value in 
Tablel^J. This makes the distance smaller by ~4 per cent. On 
the other hand, correcting the observed ffuxes for the fact 
that the RXTE/PCA ffuxes are generally found to be 20 
per cent higher than fluxes determined using other satellites 
would make the radius expansion burst distance larger. Fur- 
thermore, if Cyg X-2 is a halo object with a low metallicity 
the absolute V-band magnitude would be smaller than that 
of a star with solar metallicity in order to explain the ob- 
served spectral type, making the discrepancy even bigger. A 
possible solution for the distance discrepancy could be that 
the neutron star in Cyg X-2 is lighter than the canonical 
value of 1.4 Mq (although this would also affect the opti- 
cally determin ed distance). Howeve r, th is is at odds with 
the fin dings of ICasares et al.! il998l) and lOrosz fc Kuulkerj 
il999t) who find that the mass of the neutron star in Cyg X-2 
is >1.88 Mq and 1.78±0.23 Mq, respectively. 

From spectro scopic obs e rvatio ns of XTE J2123-058 
in quiescence jCasares et alJ (|20o3) determined a best-fit 
spectral type for the companion star of K7V (they ruled 
out spectral types e arlier than K3 and later than Ml). 
IShahbaz et al.l ll2003l) report a quiescent V-band magnitude 
for XTE J2123-058 of 22.65±0.06, a reddening in the V- 
band of 0.37 magnitudes, and that the companion star con- 
tributes approximately 77 per cent of the fiux in the R- 
band. Here we assume that the companion star contributes 
70 per ce nt of th e ffux in the V-band. Again using the re- 
lation of iPoppeil lll980i) to estimate the absolute V-band 
magnitude gives Mv = 7.97 for the observed mass ratio of 
0.49 and assuming a neutron star mass of 1.4 M0 (had we 
assumed a neutron star mass of 2.0 Mq we would have got- 
ten Mv = 7.72). We again assumed that the companion star 
fills its Roche lobe. From this we derive a distance of 8.7/9.6 
kpc for a neutron star ma ss of 1.4/2.0 Mq, respectively (see 
also iTomsick et alJ l2004l and references therein) . The dis- 
tance derived using the radius expansion burst is 15.7-21 
kpc (see Table |21 deadtime effects are negligible in case of 
XTE J2123-058). For the limiting spectral types of the com- 
panion star (K3/M0V) the distance would be 14/7 kpc for 
a 1.4 Mq neutron star. So, even for an assumed K3 spec- 
tral type the distance is lower than that derived using the 
radius expansion burst. Again, as in the case of Cyg X-2, 
correcting the parameters we used in the distance calcula- 
tions either for the fact that the RXTE/PCA ffuxes are gen- 
erally found to be 20 per cent higher than ffuxes determined 
using other satellites or for the fact that XTE J2123-058 
might be a halo object with a low metallicity would make 
the discrepancy between the photospheric radius expansion 
burst-distance and the method B bigger. 

Even though the sample of neutron star sources for 
which we can compare the distances derived using radius 
expansion bursts and that derived using the properties of 
the companion star is small (so far this is only possible for 
two source), this could mean that distances derived using 
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Table 2. Properties of the sample of neutron stars used in this Paper. Distances derived from type I photospheric radius expansion bursts 
in both persistent and transient neutron star systems (excluding those in Globular Clusters) using a neutron star Eddington luminosity 
of 2.0 or 3. 8x10"^* erg s~^. The burst peak flux (unabsorbed, 0.1—100 keV; indicated with a i behind the flux) or the bolometric peak 
burst flux (indicated with a ii behind the flux) is given as well. The z— values are rounded to the nearest 25 pc. 



Name 


1 
1 


u 
D 


1/p 


^ orb 

hours 


Peak burst flux 

— 2 — 1 

erg cm s 


Distance 

Z.U— O.O t^KpCj 


2.0-3 


z 

8^= (pc) 


References 


F,XO 074S— (S7fi 


279.98 


-19.81 


T 


3.82 




6.8-9.1 


-2300 


- -3075 


L-'-l 


2S nq 18-54 


275.85 


-3.84 


p 


7 


9.4xl0~* (ii) 


4.3-5.8 


-300 


- -400 


[18,19] 


Cir X-1 


322.12 


0.04 


T /P? 


398 


1 9x1 l{\ 


7.8—10.5 





- 




41 1 1 finS— 522 


330 93 


-0.85 


T 


12.9? 


1 2xin~^ l{\ 


2.8-3.8 


-50 


- -50 


141 


Sco X-1 


350 09 


23.78 


p 


18.9 


parallax 


2.8+0.3 


1125 


[201 


4TI 1fi'?fi-5'? 


332.91 


-4.82 


p 


3.80 


1 '^y^o~^ ({) 


3.7-4.9 


-300 


- -400 


1211 


41 J 1658-2Q8 


353 83 


7.27 


T 


7.11 


2.5xl0~* (ii) 


8.4-11.3 


1075 


- 1425 


151 


4U 1702-429 


343.89 


-1.32 


p 


7 


6.6x10"*''' (ii) 


5.3-7.1 


-125 


- -175 


1321 


4U 1705-44 


343.32 


-2.34 


p 


7 


3 7x10"*''' (ii') 


7.2-9.6 


-300 


- -400 


1321 


XTE J1710-281 


356.36 


6.92 


T 


7 


8.6x10"^'*' (ii) 


14.8-19.8 


1800 


- 2375 


[321 


SAX J1712. 6-3739 


348.93 


0.93 


T 


7 


5.1x10^® (ii) 


5.9-7.9 


100 


- 125 


[61 


IH 1715-321 


354.13 


3.06 


P/T7 


7 


6.7x10^® (ii) 


5.1-6.9 


275 


- 375 


[221 


RX J1718. 4-4029 


347.28 


-1.65 


p/T? 


7 


4.3xl0~* (i) 


6.4-8.6 


-200 


- -250 


[231 


4U 1728-34 


354.30 


-0.15 


P 


7 


8.4xl0~* (i) 


4.5-6.1 


- 


--25 


[24,25] 


KS 1731-260 


1.07 


3.65 


T 


7 


6.3xl0~* (ii) 


5.3-7.1 


325 


- 450 


[71 
I'J 


4U 173.5-44 


346.05 


-6.99 


p 


7 


2.9x10"*'*' (ii) 


8.0-10.8 


-975 - 


--1325 


[321 


GRS 1741.9-2853 


359.96 


0.13 


T 


7 


4.0x10"* (i) 


6.6-8.9 


25 


- 25 


[8] 

L J 


2E 1742.9-2929 


359.56 


-0.39 


T/P? 


7 


3.7x10"*''' (ii) 


6.9-9.2 


-50 


- -75 


[32] 


SAX J 1747. 0-2853 


0.21 


-0.24 


T 


7 


3.2x10"* (ii) 


7.5-10 


-25 


- -50 


[9] 


GX 3+1 


2.29 


0.79 


P 


7 


9.3x10-* (i) 


4.3-5.8 


50 


- 75 


[26] 


SAX J1750.8-2900 


0.45 


-0.95 


T 


7 


6.4x10"* (ii) 


5.2-7.0 


-75 - 


- -125 


[10] 


SAX J1752.3-3138 


358.44 


-2.64 


P/T? 


7 


2.8xl0"*(ii) 


7.9-10.6 


-375 


- -475 


[27] 


SAX J1808.4-3658 


355.38 


-8.15 


T 


2.0 


2.5x10-^ (ii) 


2.7-3.6 


-375 


- -500 


[11] 


SAX J1810.8-2609 


5.20 


-3.43 


T 


7 


7.0x10-* (ii) 


5.1-6.8 


-300 


- -400 


[12] 


4U 1812-12 


18.06 


2.38 


P 


7 


1.5x10"^" (ii) 


3.4-4.6 


150 


- 200 


[28] 


XTE J1814-338 


358.75 


-7.59 


T 


4.27 


2.6x10"* (ii) 


8.2-11.0 


-1075 


- -1450 


[13] 


GX 17+2 


16.43 


1.28 


P 


7 


1.2x10"* (ii) 


11.9-16.0 


275 


- 350 


[29] 


Ser X-1 


36.12 


4.84 


P 


7 


2.1x10"*''' (ii) 


9.5-12.7 


800 - 


- 1075 


[32] 


Aql X-1 


35.72 


-4.14 


T 


19.0 


9.1x10"*''' (ii) 


4.4-5.9 


-325 


- -425 


[14] 


4U 1857+01 


35.02 


-3.71 


T 


7 


3x10"* (ii) 


7.5-10 


-500 


- -650 


[15] 


4U 1916-053 


31.36 


-8.46 


P 


0.83 


3.1x10"* (ii) 


7.5-10.1 


-1100 


- -1475 


[30] 


XTE J2123-058 


46.48 


-36.20 


T 


5.96 


7x10"^ (ii) 


15.7-21 


-9275 - 


-12.4x10^ 


[16,17] 


Cyg X-2 


87.33 


-11.32 


P 


236.2 


1.35xlO"*(ii) 


11.4-15.3 


-2250 


- -3000 


[31] 



T= transient, P=persistent 

'' Corrected for the fact that RXTE fluxes are found to be higher by 20 per cent 

-1 



Determined assuming an Eddington peak flux of 2.0 or 3.8x10^* 



erg s 



Referenc es: [lllGottwald et al.lll986l [2llTennant et al'!||l986l [3llBra.ndt et al.lll99fiL [4llMurakami et alJll987l [51 IWiiiiands et al.ll200l 
[6llCocchi et alJl20 01a. [71 Muno et al.' 2001. [81 "Cocchi et al."l999'. [91 'Nat alucci et alj|2000 b'. flOl Kaaret ct al. 2002, [11] 
lin't Zand et al.|200]1 7l2l Natalucci et aL 200Qa. [1 31 .Strohmavcr et al. ,200a. [14] This work [15] X^hcvalier & Ilovaiskv 1990a, [16] 
l.lll99gl. [17llTomsick et alj|l999ar[18l Ijonker et alj|200ll [191 ICornelisse et al.ll200l 



iHoman et al.|[l99a. [17 1 iTomsick et'ar]ll999al." '[18l Ijonker et alj 200 ll ?19ll 
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super 



[20] Bradshaw ct al. 1999, [21] 
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burst" , [251 iHoffman et" al. 1980, [ 261 iKuulkers fc vaii der Kli j|20od. [271 ICocchi et alJl2001bl. [281 ICocchi et al.ll20ol [29] 



lallowav et al. II2OOII [311 ISmalelll998l ~i32l Galloway et al. 2004 
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method B are too low, or that the distance derived using 
the radius expansion burst is too large. Below we will in- 
vestigate in some detail possible effects responsible for the 
observed discrepancy in distance estimate between method 
B and the radius expansion burst method. 

2.3.1 Systematics in the companion star radius 
determination and residual accretion 

As mentioned earlier, under the assumption that the 
companion star fills its Roche lobe a first estimate 



of the radius of the compan ion star can be obtained 
(cf. IWhite fc van Paradiislll99(t) . Given the fact that mass 
accretion must take place in order to explain the multiple 
outbursts for some of the systems the assumption that the 
secondary (nearly) fills its Roche lobe seems justified. How- 
ever, if the temperature distribution on the surface of the 
star is uneven e.g. due to effects of irradiation, then the lu- 
minosity is not determined by the full size of the companion 
star. Depending on whether the determined temperature is 
that of the hotter or colder part of the star this yields an 
over or under estimate of the distance. 
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An under estimate of the amount of light contributed 
due to residual accretion would lead to an under estimate of 
the distance and vice versa. 



2.3.2 Systematica in the companion star temperature 
determination 

Besides the radius, the temperature of the companion star 
is important for its luminosity. The temperature is de- 
rived from the observed spectral type. A systematic mis- 
classification of one spectral type for a fixed radius (e.g. us- 
ing a KIV instead of a KOV star) already results in a dis- 
tance error of ~15-25 per cent for late type stars (see Fig- 
ure0. The difference in derived distance between using the 
surface brightness of main-sequence stars or giants (for fixed 
radius) is also given in Figure The giant surface bright- 
nesses lead to smaller absolute magnitudes, i.e. smaller dis- 
tances. 

The best-fit spectral type can be determined using the 
optimal subtraction method iMarsh et alj Il994) which is 
quite robust. However, this method is rarely applied fully. 
The broadband spectral energy distr ibution can also be used 
to determine the spectr al type (cf. iGelino et al.ll200ll and 
iGelino fc Harrisonl20ol for recent use of this method) . How- 
ever, disentangling the effects of fast rotation, reddening, a 
possible accretion disc contribution, and the spectral type 
is difficult with only a few data points. This makes the un- 
certainty in the temperature of the companion star an im- 
portant contributor to the uncertainty in the distance. For 
XTE J2123-058 a spectral type of K2 (just outside the deter- 
mined range) would already make the two different distance 
estimates consistent. For Cyg X-2 a spectral type of AO is 
needed while the determined type is A9±2. However, from 
the effects discussed above, there is no reason to expect a 
systematic underestimation of the effective temperature of 
the companion star which would lead to a systematic under- 
estimation of the distance using method B. 

In using method B as described above an implicit as- 
sumption about the temperature of the companion star has 
been made. This is due to the fact that the effective temper- 
ature of the companion star is in most cases assumed to be 
that of a main-sequence star of the observed spectral type. 
Besides a small variation in temperature with luminosity 
class there is also a difference in bolometric correction for 
stars which have the same spectral type but a different lu- 
minosity class. Comparing the observed surface gravity with 
that of different luminosity classes of the observed spectral 
type the effective temperature of the companion star can 
be determined. These then determine the bolometric cor- 
rection, and thus the distance. However, the error in the 
distance from neglecting the difference in bolometric correc- 
tion between the luminosity classes is small. 

However, it is possible that the observed spectral types 
are systematically shifted to later spectral types because of 
effects on the lines and continuum induced by the fast ro- 
tation of the companion star. Single late type stars have 
rotational velocities of less than a few kilometres per second 
due to the onset of a magnetic-brake in stars later than typ- 
ically F5 (see [Kraft 196^. However, the late type compan- 
ion stars in LMXBs are likely in co-rotation with the orbit. 
Therefore, the rotational velocities range from several tens of 
km s~^ for long orbital period systems (e.g. GRS 1915-f 105; 
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Spectral type 

Figure 1. The fractional error one would make in the distance 
estimate if the spectral classification is off by 1 sub-class (grey 
squares) is plotted as a function of spectral type for late type 
stars. The fractional error in the distance estimate if one takes 
a main-sequence spectral class companion instead of a giant of 
the same spectral type for the same radius (black, filled squares). 
For both it is assumed that the radius of the companion star is 
known and accurate (in the calculations we used the results from 
Popper 1980). 



iHarlaftis fc Greineilliooil to hundred or more km s ^ for 
the short orbital period systems. Studies of rapidly rotat- 
ing early type stars found that a fast rotation causes a 
slight increase in absolute magnitude of stars with spectral 
type later than B5 (the stars are intrinsically less luminous 
than the non-rotating stars of the same spectral type by 
typically several tenths of a magnitude; only stars up F8 
were studied; CgUms & Sonncbornl ll977r) . Additionally, it is 
known that limb and gravity darkening effects on the line 
change the equi valent widths of stella r photospheric absorp- 
tion lines (e.g. IShain fc Struvelfl92gl : ICollins fc SonnebornI 
Il977l) .^ Generally, the lines in the spectrum resemble the 
lines of a later spectral type star (Collins & SonnebornI 
^23) • However, this depends on the behaviour of the equiv- 
alent width of the used lines as a function of the effective 
temperature. Broadening of template star spectra in the op- 
timal subtraction technique does not take into account these 



^ e.g. the ratio between equivalent widths of the hydrogen lines 
and those of weak metallic atoms also changes among other things 
due to the fact that the strong hydrogen lines respond d i fferent 
to the rotational effects ICoUins et alJll99lft . lOrosz et all J200j) 
found that the companion star in V4641 Sgr (a.k.a. SAX J1819.3— 
2525) is rapidly rotating (v sini ~123 km s^^). In fact, for an 
assumed B9III spectral type observed at an inclination close to 
60° tOrosz et al.. .20011) . t he star is r otating at 80-9 per cent 
of its break-up speed fcf. ICoUins fc Sonncborni ll977tl . The fact 
that the lines of stars with large rotational velocities are similar 
to th e lines in stars of later spectral types <Collins fc SonnebornI 
[l977|) can explain th e unusu ally strong Mg II line in V4641 Sgr. 
iBildsten fc R,utledgel j200Cl) noted that the peculiar Li abun- 
dances observed in the companion stars of several BH SXTs could 
also result from the fact that the companion stars are fast rotators 
similarly to RS CVn stars. 
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physical changes in the hnes. In early type stars the decrease 
in luminosity of the star due to rotation is outweighed by 
the apparent shift in spectral type. Hence, the distance es- 
timation method B using spectra of non-rotating template 
stars to estimate the spectral type applied to fast rotating 
early type stars would likely lead to a net underestimation 
of the distance, the effect being bigger the later the spectral 
type of the star (at least up to spectral type F8). 

Unfortunately, there are no detailed model atmosphere 
calculations showing the physical effects of rotation on the 
absolute magnitude and the lines of stars later than F8 avail- 
able in the literature. However, effects of limb darkening are 
larger in late type stars than in early type stars. Further- 
more, extrapolating the findings for the early type stars the 
spectral appearance is likely more affected than the true 
source luminosity meaning that the distances to the sources 
would be underestimated using non-rotating template stars 
to determine the spectral type. 



2.3.3 Systematics in the interstellar absorption estimates 

If the interstellar extinction, Ay, is systematically overesti- 
mated this would yield an underestimate of the distance 
or vice versa since mv — Av — Mv ~ 51ogd — 5. The in- 
terstellar extinction caused by scattering of light off dust 
grains is traced by the optical/UV spectrum (through 
E(B-V) and/or the strength of the absorption feature at 
2175A). Spectral fits to the X-ray data trace interstel- 
lar absorption caused by neutral hydrogen and absorp- 
tion (above ^0.5 keV) by K/L-shel l electrons of mostly 
O and Fe JPredehl fc Schmitd ITqQsIi . Finally, intersteUar 
absorption lines from ions such as Na I and Ca II trac e 
low-ionisation interstellar gas dSembach fc Dankj Il994l . 
Over the years relations between Ay and these differ- 
ent sources of interstellar extin ction have been established. 
E.g. iPredehl fc Schmitti (ll995|) r eporte d a relation between 
Nh and Av. iRieke fc Lebofskl give a relation be- 

tween E(B - V) and Av (R x E(B - V) = Av, with Rr^S.I). 
[HcrbiE (1995) gives relations between the equivalent widths 
of the Na I and Ca II absorption lines and E(B — V). 
However, the line-depth of some absorption lines satu- 
rates quickly leading to an overestimate of the distance. 
Furthermore, for systems out of the Galactic plane the 
equivalent width of the interstellar absorption features does 
not increase much with distance once the system has a 
z-distance larger than the scale-height of the interstellar 
material (e.g. the scale-height of Ca II is approximately 
800 pc whereas that of Na I is approximately 400-500 pc; 
ISembach fc DanksllTggi) . 

In addition, determining Av from the observed Nh or 
from the observed equivalent width of interstellar absorp- 
tion lines one implicitly assumes that the sight-line for the 
SXT under investigation has the same gas-to-dust ratio 
(i.e. Nh/E(B — V)) as the sight-lines for which the rela- 
tions between Ay and Nh or between the equivalent widths 
of the Na I and Ca II absorption lines and E(B — V) have 
been established. It is known that this is not always the case; 
e.g. the sight-line t owards the Gould belt (i .e. Orion, Ophi- 
uchus, and Perseus; iBurstein fc Heile3ll978l) has a low gas- 
to-dust ratio. Similarly, the sight-line towards the Galactic 
Centre, where most LMXBs are located, has, in general, a 



low R ('R<3.1: IUdalskill2003l) . This caveat has to be kept in 
mind. 

In Table 13 we show the Ay as derived fro m optical/UV 
and X-ray data. As was found before (e.g. IVrtilek et alJ 
fl99lll . the extinction derived using X-ray data is system- 
atically larger than that based on optical/UV data. There 
are two possible reasons for this. iPredehl fc Schmitd 
used the Nh values derived from fitting an absorbed power- 
law model to ROSAT X-ray data. Using an absorbed 
Bremsstrahlung or a black body model to fit the data gave 
systematically lower values for Nh. However, in the Nh esti- 
mates in Table|3in most cases a power-law model has been 
used as well. 

Another possibility is that local absorbing material is 
present during outbursts. Since most of the X-ray observa- 
tions which are sensitive enough to determine Nh are done 
during outburst this could lead to a systematic overestima- 
tion of Nh. The fact that for XTE J1859+226 the value 
derived from X-ray spectral fits is higher than that derived 
using the UV data even though both the HST UV and the 
BeppoSAX X-ray observations were made during outburst 
can be explained by the fact that the UV and X-ray observa- 
tions are sensitive to different sources of interstellar extinc- 
tion. If the extra absorbing material is indeed local to the 
source, the temperature may be too high for dust to form, 
on the other hand neutral hydrogen may well be present ex- 
plaining the difference in derived Ay. This all suggests that 
the Ay derived from Nh from X-ray spectral fitting, sys- 
tematically provides an overestimate for Ay. Hence, their 
use would yield an underestimation of the distance. 

In all BH distance estimates in the literature the Ay de- 
rived using optical or UV data has been used, except in case 
of XTE J1550-564, H 1705-25, an d GS 2000-K25. The diffuse 
interstellar bands (DIBs) used bv ISanchez- Fernandez et alJ 
|1999) to estimate the Ay for XTE J1550-564 do not suf- 
fer from saturation at relatively low values of E(B — V). 
Given the fact that the factor with which the distance to 
the source changes goes as ^ = 10 » , the distance 
to XTE J1550-564 may have been underestimated by a fac- 
tor ~3! Although, as explained above, the sight-line could 
have properties different from the properties of the sight- 
line used in deriving the relation between the Ay and the 
equivalent width of the DIB leading to an anomalous low Ay 
as derived from the DIB. In case of GS 2000-1-25 an Ay of 
4.4 has been used. However, since the optically derived value 
for Ay for GS 2000-f 25 does not seem to be very accurate 
(it has been determined assuming the intrinsic B-V colour 
of the source during outburst is 0) it is unclear whether this 
is an overestimation or not. Unfortunately, it is not possible 
to estimate whether the distance to H 1705-225 is over or 
underestimated since there is no optically derived Ay avail- 
able. 



2.3.4 Systematics in the radius expansion burst method 

It has been suggested that the burst flux could be anisotropic 
|Kuulkcrs et al. 200^. However, given the good agreement 
between the Globular Cluster distances and the distances to 
the LMXBs in these Globular Clusters derive d from the pho- 
tosph eric radius expansion burst properties ^^jffke^^^I 
l2003r that seems unlikely. Furthermore, ICallowav et al. 
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Table 3. The interstellar extinction Ay for the sample of BH 
SXTs determined using optical/UV or X-ray data. For some 
sources a range of values corresponding to the different values 
found by the different references is given whereas for others an er- 
ror bar is given. We used R = 3.1 in Ay = R X E(B — V). We used 
the Njj value obtained using a power— law model where available 
since the correlation Nfj = 1-79 X 10'^^ cm~^ Av was derived us- 
ing power-law fits to ROSAT data (see Predehl & Schmidt 1995). 
An "X" denotes that the value has not been determined. 



Name 


Av optical 


Av X-ray 


Rcfs 








Opt - 


X-ray 


GRO J0422+32 


0.6-1.2" 


<2.8 


[1]- 


[2] 


lA 0620-00 


1.09-1.24 


1.3±0.7 


[3,4] 


-[5] 


GS 1009-45 


0.6±0.2 


X 


[6]- 


[X] 


XTE J1118+480 


X 


0.041^10.004"^ 


[X]- 


-m 


GS 1124-684 


0.9±0.1 


1.28±0.06 


[8]- 


[9] 


4U 1543-47 


1.55±0.15 


2.4±0.1 


[10]- 


[11] 


XTE J1550-564 


2.5*= 


4.88±1.15 


[12]- 


[13] 


GRO J1655-40 


3.7±0.3 


4.8±2.8 


[14] 


-[5] 


GX 339-4 


>2.8 


3.9±0.5 


[15]- 


[16] 


H 1705-250 


X 


1.7±0.5 


[X]- 


[17] 


SAX J1819.3-2525 


1.0±0.3'' 


X 


[18]- 


-[X] 


XTE J1859-f226 


1.80±0.37 


4.47*= 


[19]- 


[20] 


GS 2000-f 25 


3.5'' 


6.4±1.0 


[21]- 


[22] 


GS 2023-f338 


3.3-4.0 


3.9±0.4 


[23,24 


-[5] 



° see the discussion and references in lBeekman et all il997ft 
the uncertainty is large since the value is derived assuming 

{B - y)o = 
no error bars given 
Nh from EUVE observat i ons 

" ISanchez-Fernandez et alJ ^1999^ report 2.2±0.3. We prefer the 

value derived using the DIB since the Na D line may have been 

saturated. 



Re ferences: fl] iBeekman et al.l Il997t [2] [Sunvaev et al.l Il993t 



[31 lOke fc GreensteinI 11977 
20021: [61 Idella Valle et al " 



[41 IWu erai]_ 19831 [5] [ Kong et aO 
119971: [7] iHvnes et al.1 l200d 



Cheng et al][l9 92': [9] 'Grciner et al."l994'; [10] "Orosz et al."l998'; 
]ll]_ lvan der Wocrd ct al. 1989; [12] Sanchez-Fernandez ct al. 
1999; [13] Tomsick ct al. 2001; [14] Hvncs ct al. 1998; [15] Hynes 
et al. (2004); [16] .Gallo ct al...200a: [17] .Griffiths ct al.a973 [181 



Orosz et al.l l200li: [19] iHvnes et 311120021: [20] Idal Fiume^t al.1 
19991 [21] IChevalier fc Ilovaiskvlll990d: [22] rTsunemi et al.lll989t 
[23] IWagner et alj|l99ll: [24] ICasares et aljil993l 



i2003D showed that considering the burst peak fluxes of 
61 photospheric radius expansion bursts in the atoll source 
4U 1728-34 the degree of anisotropy in the burst emission 
in less than 2 per cent. 

Two out of the eight neutron star systems studied by 
iKuulkers et al . (2003) have a photospheric radius expansion 
peak burst luminosity that is lower than that of the other six. 
This lower peak luminosity is consistent with the Eddington 
luminosity limit for hydrogen rather then helium-rich ma- 
terial for a neutron star mass of 1.4 IVIq. So, for some of the 
neutron star systems in our list we could have overestimated 
the dista nce by a factor ^ Vl. 8. If we take a 1:4 ratio as was 
found bv lKuulkers et alJl200.'j this would affect 8 sy stems in 
our sample. As was noted bv lKuulkers et al.l ll2003l) the two 
sources whose peak luminosity is consistent with hydrogen 



accretion have an orbital period characteristic for normal 
LlVEXBs, rather than ultra-compact X-ray binaries with pe- 
riods less than ~1 h our (which n ecessarily accrete hydro- 
gen poor material e.g. lVerbunt fc va n don Houvel 1995). Be- 
cause the fraction of ultra-compact X-ray binaries in Glob- 
ular Clusters might be higher than in the field, the ratio 
between the number of sources with a hydrogen-rich, low 
Eddington burst luminosity and the number of sources with 
a helium-rich, high Eddington burst luminosity could be 
higher. It is unlikely, however, that the distances to all neu- 
tron star LlVIXBs in our sample have been overestimated. 
E.g. in the long period LMXBs KS 1731-26 and 4U 1636-53 
the photospheric radius e xpansion burst was consistent with 
a hel ium-rich explosion i]\fa^i£_et_jjJ|200CjJSu ^ all 
Il984l respectively, see also lCumming fc Bildstenll2000l) . 



3 IMPLICATIONS AND DISCUSSION 

Using the data set compiled above, leaving the distances 
derived using method B as they are, keeping the discrepancy 
in the distance estimate between method B and the radius 
expansion burst method in mind, we investigate the Galactic 
distribution of neutron star and BH LMXBs and the peak 
luminosity of BHs. 



3.1 Galactic distribution of LMXBs 



Ivan Paradiis fc White! lll995l) investigated the Galactic z- 
distribution of neutron star LlVEXBs and concluded that neu- 
tron stars should receive an asymmetric kick at birth from 
the fact that the rms value of the z-distribution was ^ 1 
kpc. We obtain an rms z-values of 1025 pc and 1125 pc for 
persistent and transient neutron star LMXBs, respectively 
(we round rms z-values to the nearest 25 pc). In this we have 
excluded XTE J2123-058 since with a z-value of -1 2 .4 kpc 
it would dominate the outcome and ICasares et al.l (|2003) 
show that the systemic velocity is consistent with it being a 
halo source. It can be argued on the basis of the large z val- 
ues that Cyg X-2 and EXO 0748-67 6 are also halo sou rces 
although this is much less clear. E.g. iKolb et alj ll2000l) ar- 
gue that if Cyg X-2 is at a distance of ~11.6 kpc it could 
have originated in the Galactic plane. However, when we 
also exclude Cyg X-2 and EXO 0748-676 we find an rms z- 
value for the persistent and transien t neutron star LMXBs 
of 700 pc and 850 pc, respectively. Ivan Paradiis fc Whitd 
il995h found an rms-z value of ~500 pc when they excluded 
Cyg X-2 and EXO 0748-676 from their sample. Since we 
have used 3.8x 10^* erg s"'^ for the Eddington luminosity for 
the radius expansion peak luminosity for all neutron stars, 
the rms z-value we derived for neutron star LMXBs cor- 
responds to an upper limit. If we use 2.0x10^* erg s~^ as 
the photospheric radius expansion peak burst luminosity we 
find an rms z-value of 775 pc and 850 pc for persistent and 
transient neutron star LMXBs, respectively (here we only 
excluded XTE J2123-085, if we also exclude Cyg X-2 and 
EXO 0748-676 we find 550 and 650 pc for persistent and 
transient neutron star L MXBs, respectively). I. e . our find- 
ings confirm the result of Ivan Paradiis fc WhiteHl995l) 

In a follow-up paper, IWhite fc van Paradiij ITggg) m- 
vestigated the differences between the rms values of the 
neutron star and BH LMXBs z-distributions. They found 



that the rms value for the BHs was substantially lower than 
that of the neutron star z-distribution (more than a fac- 
tor 2). Using the distances of the BHs given in Table Q 
we now find an rms-value of ~625 pc (we took a distance 
of 2.5 kpc for GRO J04 22+32; we excluded the hkely halo 
object XTE ,71118-1-480: IWagner etalltoll [ahhough note 
that some of the evidence leading to the suggestion that 
XTE J 1118-I-480 is a ha lo object was based on the low rms 
z-value lwiiite fc van Paradiis 199ft found for BH LMXBs]), 
i.e. close to the upper limit we find for the neutron star 
systems. Increasing the distance of XTE J 1550-564 with a 
factor 3 does not significantly increase the rms value of the 
BH population. The main reason for the diff erence between 
the findings of IWhite fc van Paradiij il996l) and our find- 
ings is that the distance estin iates for most BHs have gone 
up. The conclusion drawn bv lWhite fc van Paradii^ lll996ll 
that BHs receive a significantly smaller kick than neutron 
stars is no longer tenable. 

We plotted the z values for transient and persistent neu- 
tron star LMXBs (open diamonds and squares, respectively) 
and BHs (black dots) using the distances from Tableland 
|2| as a function of their projected distance to the Galactic 
Centre in Figure |21 It is interesting that the rms value for 
the projected distance to the Galactic Centre for the neutron 
stars and BH is 4.8 kpc and 7.0 kpc, respectively (excluding 
the [likely] halo sources). Because the Galactic potential in 
the z-dir ection decreases with increa sing Galactocentric ra- 
dius (e.g. ICarlberg fc Innanenlll98'^ the larger scale height 
of BHs can partly be due to this eff ect, r ather tha n a kick 
veloci ty. E.g. for neutron star svstems lvan Paradiis fc Whit j 
(^2^ estimate an rms value of the z-distribution near the 
Sun of 650 pc. Furthermore, the symmetric kick velocity 
(imparted due to mass loss in the supernova) scales with 
the companion mass and the mass lost in the supernova 
(e.g. Nclcmans et al. 1999) both can be (much) larger in the 
case of BH systems, but it scales inversely with the total 
mass of the remaining binary. Lastly, for large asymmetric 
kick velocities the binary is more likely disrupted in case of 
a neutron star than a BH. A detailed investigation of all the 
possible explanations for the high rms value for the BHs is 
needed in order to draw firm conclusions and we defer this 
to a later paper. 

After lwhite fc van Paradii j l)l996h we also compare the 
distribution of the 1 and b coordinates of the neutron star 
and BH LMXBs. We first compare the 1 and b coordinates 
of neutron star LMXBs in Table |21 with those of the BHs 
listed in Table A K-S test shows that the probability 
that the distributions are the same is 18 per cent for the 
I coordinate and 52 per cent for the b coordinate (the K- 
S D-value is 0.33 and 0.24, respectively. Here and below: 
the effective number of data point is always larger than 
10 so that the probabilities we quote are quite accurate, 
cf.^rcss et al. 1992j. To investigate this further and to min- 
imise selection effects, we plot the 1 and b coordinates of all 
neutron star systems , e.g. systems wh ere a burst was found, 
and pulsars listed in iLiu et al.l i200ll) (solid line histogram 
in Figure ISJ. We compare this with the BH sources listed 
in Table Q and the BH candidates for which no mass es- 
timate bas ed on a radial velocity study exist (systems in 
table 3 in iMcChntock fc R.emiUardI l2004l : the dashed line 
histogram in Figure O. A one-dimensional K-S test shows 
that the probability that the neutron star and BH distri- 
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Figure 2. The z-distribution of the Galactic persistent (open 
squares) and transient (open diamonds) neutron star LMXBs 
for which a radius expansion burst has been detected, and BH 
LMXBs (filled circles) for which a dynamical mass has been de- 
rived as a function of the projected distance to the Galactic Cen- 
tre. The projected distance is y'lP + y^ , where x and y are the 
Cartesian coordinates in the Galactic plane (see Figure The 
location of the Sun at an assumed distance of 8.5 kpc is indicated 
with an open circle, and the location of the Galactic Centre is 
indicated with a cross. Four neutron star systems fall off the Fig- 
ure (EXO 0748-676, z=-3.1 kpc, XTE J2123-058, z=-12.4 kpc, 
XTE J1710-281, Dec =19.7 kpc, and Cyg X-2, Dqc =15 kpc). 
In this Figure we use the neutron star distances derived assuming 
the Eddington peak burst luminosity was 3.8x10^* erg s~^. 

butions in 1 are drawn from the same distribution is 37 
per cent whereas the probability that b values are drawn 
from the same distribution is 90 per cent (the K-S D-value 
is 0.19 and 0.12, respectively). Hence, using a larger sam- 
ple of neutron stars and BHs makes the probability that 
the neutron star and BH distributions are the same larger. 
Th e probabilities we derive are higher than those derived 
bv lWhite fc van Paradiis' (1996'); a difference we attribute 
to the fact that White fc van Paradiis (1996) had fewer BH 
systems in their study. We conclude that there is no evidence 
for a difference in the 1 and b distributions of neutron star 
and BH LMXBs. 

Interestingly, there seems to be an excess of 12 out of 
a total of 22 BH and an excess of 19 out of a total of 45 
neutron star systems with -30° < 1 < 0° over systems with 
0° < 1 < 30°. There is an excess of 5 out of 35 BH and 
13 out of 53 neutron star systems with -10° < b < 0° and 
0° < b < 10° (see Figure EJ. The expected mean number 
of systems with a positive/negative 1 of a symmetric dis- 
tribution around 1=0° is N/2, with A'^ the total number of 
systems. The variance on this number is N/A. Comparing 
this with the observed number of systems yields significances 
of the observed asymmetries of 2.8cr for the neutron star dis- 
tribution in 1, 3.0(7 for the BH distribution in 1, and 3.8(t for 
the combined neutron star and BH distribution in 1. The 
combined significance for the neutron star and BH distribu- 
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Figure 3. Left panel: The 1 distribution of the neutron stars (solid line histogram; all sources showing bursts or pulsations listed in 
Table|51and Liu et al. 2001) in bins 10° wide. The 1 distribution of the BHs listed in TableQ]and McClintock and Remillard 2004 (dotted 
line histogram). Right panel: The b distribution for the same sources and using the same symbols as in the left panel in bins 5° wide. 



tion in b is l.Qcr only, so we do not discuss this apparent 
asymmetry in b further. 

The asymmetry in the 1-distribution is not an obvious 
observational selection effect. For instance, the BeppoSAX 
WFC monitoring the Galactic Centre region had a field-of- 
view of 40° X 40° centred on the Galactic Centre l|in't Zandl 
I2OOII) . The dust maps from ISchleeel et al.l Jl998l) do not 
show a large dust asymmetry around the Galactic Centre. It 
seems that the neutron star and BH LMXB distribution in 
1 is symmetric around l=-5°. Possibly this asymmetry is re- 
lated to the bar structure in the inner part of our Galaxy. It 
is known that the bar causes asymmetries in the stellar and 
gas distributions around (l,b) = (0,0) (e.g. see the recent 
review bv iMerrifield 2004). Furthermore, Gilfanov ( 2004) 
found that the distribution of LMXBs follows the stellar 
mass distribution, and hence not the spiral arm structure. 

Using the distances listed in Table Q and |21 we also plot 
in Figure^Jthe Galactic x-y distribution of the neutron stars 
for which a photospheric radius expansion burst was ob- 
served and the BHs for which a dynamical mass estimate 
has been derived. The Galactic spiral arm structure accord- 
ing to the model described in iTavlor fc C ordcs (1993) has 
been overplotted. It is clear that few LMXBs at the other 
side of the Galactic Centre have been detected, especially 
if the Eddington limit for hydrogen rich material for a 1.4 
M0 neutron star is applicable to the peak burst flux of most 
photo spheric radius expansion bursts (see also Grimm ct all 
|2002). There also seems to be a pau city of nearby neutron 
star LXMBs (this was also noted bv lVerbundbOOlD . In or- 
der to try to quantify this we per formed a two dim ensional 
Kolmogorov-Smirnov (K-S) test JPress et al.lll992D to test 
whether the neutron star and BH populations are the same. 
The K-S test gives a D of 0.45 and a probability that the 
two populations are the same of 3.8 per cent. If we decrease 
the distance of all neutron star systems by a factor ~ Vl-8 
the probability that the distributions are the same increases 



only slightly to 5.2 per cent (D=0.44). It is possible that the 
discrepancy in the spatial distribution of the neutron stars 
and BHs is a selection effect. Since we only included BHs for 
which a mass limit larger than 3M0 has been derived, the 
optical counterpart must have been detected in quiescence, 
hence this favours nearby systems. 

3.2 The black hole outburst peak luminosity 

The BH distances are important for the maximum observed 
BH luminosity. We found a somew hat larger distance for 
V404 Cyg than often used previously. IShahbaz et alJ il994lj) 
found an upper limit on the distance of V404 Cyg of 3.7 kpc 
assuming the peak outburst source luminosity was Edding- 
ton limited taking their 90 per cent confidence upper limit 
to the mass of the BH of 15 M©. The distance of 4:.Qt\% 
kpc we find shows that the maximum observed luminosity 
exceeds the Eddington luminosity for a 10 M0 BH (the limit 
often quoted to decide whether a source is a ULX or not). 
The distance derived for XTE J1550-564 may have been 
underestimated by as much as a factor 3: its distance could 
be as large as 3x5.3 = 15.9 kpc, although the uncertainties 
are large. For this distance the outburst peak luminosity in 
the 2-20 keV band al one would be 7x10' ^'^ erg cm"^ s"^ 
(taking the flux from ISobczak et al.llT999l) . For 4U 1543- 
47, SAX J1819.3-2525 and GRS 1915-^105, systems for 
which the distance was determined by method B, B, and 
C, respectively, it was noticed earlier that the outburst 
peak lumi nosities are super-Eddington for the best-fit BH 
mass e s ('cf. iRevn ivtsev et al.ll20o'3 iMcClintock fc RemillardI 
I2OO4I) . [Garcia et al.l ll2003l) show that the peak outburst of 
GRO J 1009-45 most likely also exceeded the Eddington lu- 
minosity for a for ~ 10 Mq BH (even if the distance is 5.7 
instead of the 9 kpc they favoured). It seems therefore likely 
that there are several sources in our own Galaxy which we 
would classify as transient ULXs had we observed them in 
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Figure 4. Left panel: The x— y distribution of the Galactic persistent (open squares) and transient (open diamonds) neutron star LMXBs 
for which a photospheric radius expansion burst has been detected, and BH SXTs (filled circles) for which a dynamical mass has been 
derived. The location of the Sun is indicated with an open circle, and the location of the Galactic Centre at an assumed distance of 
8.5 kpc is indicated with a cross. It is clear that there is a paucity of nearby neutron star systems compared to nearby BH SXTs. The 
neutron star distances were derived assuming that the Eddington peak burst luminosity was 3.8x10^* erg s~^. Overplotted is the spiral 
structure and the free-electron density according to the model of Taylor & Cordes (1993). Right panel: same as the left panel except that 
the neutron star distances have been derived assuming that the Eddington peak burst luminosity was 2.0x10^* erg s~^. 



other Galaxies. For these sources a mass determination has 
shown that they are not intermediate mass BHs but rather 
5-15 Mq BHs. Furthermore, the fact that LMXB sources 
(an old population) seem to be capable of producing super- 
Eddington luminosities could also help explain the presence 
of ULXs in elliptical Galaxies. Assuming that most ULXs 
are stellar mass BHs the fact that many ULXs are found a 
few arc seconds away from the young stellar clusters could 
be caused by kick velocities received at BH formation. 

If super-Eddington luminos ities were to be explained 
by effects of "mild-beaming" fcf. lKing et "3111200 ll) then one 
would naively expect the inclinations of the sources with 
the highest luminosities to be lowest. In Figure |S] we plot 
the inclination distribution for 13 BH SXTs in Figure |K| 
(the inchnations for GX 339-4 and XTE J 1859+226 are not 
well known and those sources have therefore not been in- 
cluded; inclinations from ^Orpsz 2003). Although the amount 
of sources is low we see that the inclination distribution of 
the BH ULX sources is not skewed to low inclinations but 
instead clusters around the fiducial 60° point. We did not 
include XTE J1550-564 (i ~ 72 ± 5°) and GRO J0422+32 
(i-41±3°) as ULXs. 



4 CONCLUSIONS 

We studied the distance estimates for BH SXTs in detail. 
We note that the uncertainties in the distance estimates 
due to e.g. uncertainties in the spectral type (i.e. temper- 
ature) of the companion star are probably large. Comparing 
the distances derived for the neutron stars Cyg X-2 and 
XTE J2123-058 using the distance estimation method B, 
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Figure 5. Inclination distribution for the dynamically confirmed 
BH SXTs (all the inclinations are taken from Orosz 2003). The 
shaded bins indicate the five systems which would have been clas- 
sified as ULX had we observed them in another Galaxy at their 
outburst peak luminosity. XTE J1550-564 (i ~ 72 ± 5°) is possi- 
bly also an ULX (see text). 



which is used for most BH SXTs, with the photospheric ra- 
dius expansion burst method, we find that the latter gives 
larger distances. This could mean that for some reason 
method B systematically underestimates the distance. Pos- 
sibly this is related to an erroneous spectral classification 
of the companion star caused by its fast rotation. If this is 
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indeed the case this would have important consequences for 
the reported difference in quiescent X-ray luminosities of 
BH and neutron star SXTs, for the BH SXT peak outburst 
luminosities, for the BH masses, and for the Galactic dis- 
tribution of BH LMXBs. Wc find that the distance towards 
XTE J 1550-564 may have been underestimated by as much 
as a factor 3 because the interstellar extinction could have 
been overestimated in that case. As was noticed before, sev- 
eral BH SXTs observed in our Galaxy would be classified as 
ULXs had wo observed them in another Galaxy (at least 5 
but perhaps even 7 out of the 15 dynamically confirmed BH 
SXTs). This suggests that many (transient) ULXs in other 
Galaxies could well bo stellar mass BHs. A re-cvahiation of 
the distance to the Galactic plane of neutron star and BH 
LMXBs shows that there is no longer evidence for a smaller 
rms-value of the z-distribution for BH systems. Such a dif- 
ference had been interpreted as evidence for the absence 
of asymmetric kicks during BH formation. However, before 
firm conclusions can be drawn about the similarities or dif- 
ferences between neutron stars and BH kicks the details of 
the formation and Galactic distribution have to be investi- 
gated. Finally, we found that the 1 distribution of Galactic 
LMXBs is asymmetric around 1=0°. 
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